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Abstract 

A new type of a capacitor with a very high volumetric capacitance is proposed. It is based upon 
the known phenomenon of a sharp increase of the dielectric constant of the metal-insulator com- 
posite in the vicinity of the percolation threshold, but still on the insulator side. The optimization 
suggests that the metallic particles should be of nanoscale and that the distance between planar 
electrodes should be somewhat larger than the correlation length of the percolation theory and 
~ 10 to 20 times larger than the size of the particles while the area of the electrodes might be 
unlimited. The random electric field in the capacitors is found to be larger than the average field 
corresponding to the potential difference of electrodes. This random field is potentially responsi- 
ble for dielectric breakdown. The estimated breakdown voltage of the new capacitor shows that 
the stored energy density might be significantly larger than that of electrolytic capacitors while 
the volumetric capacitances might be comparable. The charging and discharging times should be 
significantly smaller than corresponding times of batteries and even electrolytic capacitors. 

PACS numbers: 71.30.+h,73.63.-b, 73.22.-f 
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Creation of a capacitor with a very large capacitance per volume or weight (volumetric 
or weight capacitance) is an extremely important step on the way to a "green energetics". 
Currently the leading role in this field belongs to the double-layer electrolytic capacitors also 
known as super- or ultra- capacitors. The main idea of the electrolytic capacitor is based 
upon a very narrow spatial gap between metallic electrode and electrolyte that plays the 
role of another electrode. The modern development, both experimental and theoretical can 
be found in the papers^- and many others. The volumetric capacitance of these devices^ 
might be about 85F/cm 3 

In this communication I propose a completely different approach to the problem of a 
supercapacitor that according to my estimates may successfully compete with the previous 
ones. The basis of the new approach is the sharp increase of the dielectric constant near 
the percolation threshold in a system of metallic nanoparticles distributed randomly in a 
dielectric matrix. 

Suppose that the permittivity of the dielectric is e K D , where e is the permittivity of 
vacuum, and is the dielectric constant. Suppose that a; is a volume fraction of metal in 
the mixture. The insulator-metal transition of the percolation nature occurs at some critical 
value x c . This critical fraction is practically independent of the form of the particles and 
x c ~ 0.15 in three dimensional systems with a random distribution of metallic particles^. 

It has been shown by Dubrov, Levinshtein and Shur 7 and by Efros and ShklovskiiS that 
on the dielectric side of the transition the dielectric constant k of the composite material 
with the infinite volume diverges near the percolation threshold as 

Kd ft \ 

K, = , (1) 

t = x c — x > 0, and q is a critical exponent. Note that Eq.flTJ) is only valid at x c — x <C x c . 
A similar result was obtained later by Bergman and Imry 9 . 

The origin of the exponent q is as follows. Consider a random composite of a dielectric 
with a small conductivity <7d and metallic particles with a high conductivity om such that 
(Jm 3> o D . One can introduce two different laws for the resulting conductivity: 



<j(x) = <7 D / (x c - X) q 



(2) 



for x < x c , and 




(3) 



2 



for x > x c . Eq. fl2]) is valid if ct.d <C cr(x) <C cr(x c ), while Eq. ([3]) is valid if <jj\/ ^> cr(x) ^> 
cr(x c ). A smooth transition from Eq.([2]) to Eq.([3]) occurs in a narrow region around x c and 8 . 

a(x c ) = a M (—) s - (4) 

It has been shown& that the three exponents q, t, s are connected by one relation 

(5) 

Eq.([T]) for static dielectric constant has been obtained by considering the same problem 
at a finite frequency u, doing analytical continuation^ in the upper half plane of a complex 
u, and then putting u — 0. The exponent q in Eq.flTJ) is the same as in Eqs. fl2fo"|) 

There are some controversies concerning the numerical value of the exponent q in three 
dimensional (3D) case. The experimental result by Grannan et al. 10 is q = 0.73. The same 
number is given in the review paper by Clerc et al.— that summarizes the results of different 
theoretical approaches. However, a similar review by Bergman and Stroud^ gives q = 0.76. 
The last two results are close and in this paper I make a choice in favor of experimental data 
and the review^, taking q = 0.73. 

The sharp increase of the dielectric constant has been confirmed experimentally in a few 
other papers^ 1 ^. In the experiment^ the dielectric constant increases about three order of 
magnitude reaching 10 5 near the percolation threshold of molybdenum particles in a ceramic. 
A similar result was reported by Pecharomman et al. 15 in BaTiC>3 — Ni composite. In the 
percolative regime they got a high and frequency independent dielectric constant k ~ 80000. 

The physical reason of such increase is the existence of metallic clusters that have large 
but finite size. This size tends to infinity near the threshold. Almost all metallic sites of 
such a cluster belong to the "dead ends" which are connected to the rest of the cluster by 
one site only. The huge capacitance is created by the small distances between sites that 
belong to different metallic clusters, separated by polarized dielectric. 

The percolation transition occurs when these clusters become connected and form the 
infinite cluster. The characteristic size of the finite clusters below the transition is of the 
order of the correlation length 

L c = 1/t", (6) 

where I is the size of a metallic particle and v is the exponent of the correlation length. In 
the 3D case u w 0.87^^. 
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Similar to Eq.flT]), Eq.(j6]) is is valid only at x c — x <C x c . Note that both these equations 
and all equations obtained from them below are rather estimates because they should contain 
unknown numerical coefficients. 

In a system of finite linear size the apparent divergence of the dielectric constant saturates 
at such small values of r, when the correlation length L c becomes of the order of this size. 
This result is common for all second-order phase transitions. 

Consider a parallel plate capacitor with a large area S of the metallic plates and with a 
distance L between the plates. The capacitance C = eqkS/L. Here k itself increases with 
L at L < L c . One can say that the dielectric constant near the percolation threshold has 
a spatial dispersion at L < L c . At L > L c it is given by Eq. (pQ) and it is independent of 
L. Then the capacitance decreases with increase of L at L ^> L c as 1/L. Thus, the optimal 
distance between the plates of a capacitor is L w L c . 

However, one should keep in mind that the percolation threshold x c has a meaning for 
the infinite array only. Any set of samples with the same size L has a small critical region 
of \x c — x\, where the threshold may appear with the probability about 1/2. The width of 
this region is&& of the order of S = (l/L) 1 ^. To be sure that there are no electrical shorts 
between metallic electrodes one should work on the dielectric side outside the critical region 
at r > 5. This means that L > L c . Since the numerical factor in Eq.Q is unknown and 
it may depend on the shape of metallic particles, I assume for the estimates that distance 
between electrodes is L c as given by Eq.([6]). To prevent electrical shorts every metallic 
electrode may be separated from infiltrated dielectric by a thin layer of a pure dielectric 
with a high electric strength. 

In this approximation the capacitance per area C$ has the form 

C s = e K/L c . (7) 

It is convenient to use a dimensionless length P such that L c — PI, where P = t~ v . Using 
Eqs. (fTTl6l) one gets 

C s = e K D (P)^- 1 ) //. (8) 
Following Clerc et al— we assume that in a 3D array the ratio qjv = 0.84. Then 

C s = e K D /(Pf 16 l. (9) 

An important characteristic of a capacitor is its volumetric capacitance. Assuming that 
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the distance between electrodes is L c one gets for capacitance per volume 

Cy = 6 K/L 2 C . (10) 



Using Eqs.( lT|6|) one gets 



or 



T \ q/u 

C v = e K D [J\ L- 2 (11) 



C v = e n D (P)(^- 2 ) jl\ (12) 



Assuming qjv = 0.84 one gets 

Cv = (13) 

The most striking result is that near the percolation threshold capacitances per area Cs 
and and per volume Cy depend on the distance between metallic electrodes much more weakly 
than in the case of a regular dielectric. Indeed, C,s ~ 1/L a16 instead of Cs ~ 1/L in a 
regular dielectric while Cy ~ 1/L 116 instead of Cy ~ l/-^ 2 , where L is the distance between 
the plates of the capacitor. This gain in the capacitance results from the spatial dispersion 
of the dielectric constant k, and it follows formally from Eq.((T]) that has been substituted 
into Eqs.f l7lfT0l) . 

Now we see that close proximity to the percolation threshold is not the optimal solution 
because in this case L c becomes larger and volumetric capacitance becomes smaller. 

One should keep in mind, however, that the strong enhancement of the dielectric constant 
near the percolation threshold is valid if L c » I which also means that P » 1 and 

T- 1 » 1. 

The computational experience shows that in a 3D array it is enough to have P ~ 15 to 
20. To find r one should use P = l/r u . 

Now I can do estimates. At P = 20 and I = 1 nm one gets 

Cs « 0.55^F/m 2 (14) 

and 



C v = 27^F/cm 3 . (15) 
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FIG. 1: Cross-section of a capacitor consisting of five thin capacitors connecting parallel. Straight 
thick lines show metallic electrodes. All space between them filled with dielectric infiltrated by 
metallic nanoparticles except thin preventing layer of a pure dielectric around each metallic elec- 
trode (see text). Thin lines show the battery connection to the capacitor. The openings between 
horizontal and vertical lines are shown large for clarity only. In fact, they should be large enough 
to prevent electrical shorts only. 

Note that the value of k^> in this estimate might be significantly larger than 100. For 
example recently™ kd = 3000 has been discovered in Sri_ x Pr x Ti03 ceramics. Taking this 
value one finds that Cy = 810F/cm 3 . Thus this estimate shows that the proposed capacitor 
might be comparable with the electrolytic one that provides Cy = 85F/cm 3 . 

The capacitor under discussion is very thin though the area of the metallic electrodes 
might be large. It is possible to make a capacitor of a cubic or of any other form with a 
parallel connection of the thin capacitors. An example of such a capacitor is shown in Fig. 
1. It has almost the same volumetric capacitance as each thin capacitor. 

Now I discuss the important issue of the maximum energy stored in the capacitor at a 
given capacitance. This energy is limited by the value of the voltage that can be safely 
applied to the capacitors without causing dielectric breakdown. 

One can show that a capacitor made from a dielectric filled by metallic particles in the 
vicinity of the percolation threshold has a microscopic field e of random orientation and 
much larger than the macroscopic field E = U/L c , calculated from the applied voltage U. 

Now I estimate the value of the field e. The energy Q of the capacitor per unit surface 
area is 

Q = C s U 2 /2. (16) 
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On the other hand, this energy can be calculated microscopically as the energy of electric 
field e per unit surface area of the capacitor 

Q = If / (G) W ' (17) 

where the integral is taken over the cube L^, but the integrand is non-zero in the dielectric 
regions only because the field inside metals is zero. Note that Eq.flTTJ) contains dielectric 
constant kd of a pure dielectric. Increase of the macroscopic dielectric constant k and of 
the capacitance is due to the microscopic field e. This field is the solution of the Poisson 
equation that takes into account metallic infiltration. Introducing the squared field averaged 
over the cube L\ 

< (e) 2 >= {l/LD J (e)W, (18) 
one can write the energy of the capacitor Q per unit area in the form 

Q = ^ < (e) 2 > L c . (19) 

Comparing Eqs. (116f8p . and (1T9]) one gets 

e « V< (e) 2 > = (L c //) ( ^" 1)/2 . (20) 

V ^c' 

Note that 

e/E « {L c /l) q/2v » 1. (21) 

Strictly speaking the field e is inhomogeneous, but the inhomogeneity is not strong. The 
microscopic field and the large value of the permittivity near the percolation threshold are 
created by the dead ends of the finite clusters. The large capacitance can be understood in 
terms of the large areas covered by the dead ends of different finite clusters with different 
potentials. This coverage is dense so that the microscopic field should not be strongly 
inhomogeneous. Therefore I think that the fluctuations of the field e are not very important, 
and this field can be used to estimate the breakdown field. 

One should also take into account that an electrical short between two finite clusters that 
appears due to a fluctuation means only that they become one cluster. Only multiple shorts 
may create the percolation between the electrodes that destroys the device. 

To estimate the breakdown voltage Uj, assume that the breakdown in the dielectric occurs 
due to tunnelling current which is usually the main reason at small distances^. The field 
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e changes its direction at a distance of a few values of / ~ lnm. Then the typical voltage 
drop inside the dielectric is of the order of 

el = l ., U . n . . (22) 

p(l-q/2u) v I 

The tunnelling current appears when the voltage drop given by the left hand side of 
Eq.( l22l) is equal to the gap width of the dielectric E g (in volts). Due to the fluctuations of 
the field e it might be somewhat less than E g , say E g /y/2. Then from Eq.( l22|) one gets 

Ub = E± P (i-«M. (23) 
V2 

At P = 20 and E g = 3V the maximum voltage that can be applied is Ub = 12V. 

These calculations reveal another important advantage of the proposed capacitor. The 
breakdown voltage is larger than E g j v2 in volts. Note that the random field e that creates 
breakdown is larger than the average field E, but the tunnelling distance I is so small that 
the voltage drop el is smaller than potential difference U that is applied to a capacitor (See 
Eq.(j22])). 

Now we can estimate the maximum density of the stored energy W as 
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W = (24) 



Using EqsflmH]) one S ets 



W = \^{E g /Wy. (25) 



This equation does not contain any critical exponents. If E g = 3V, I = lnm, one finds 
W = 0.55(/«£)/100)Wh/cm 3 . This is a very large volumetric energy. For the electrolytic 
capacitors the energy per weight is about 5 Wh/kg (See^). This happens because the 
breakdown voltage in a double-layer electrolytic capacitor is only (2 to 3)V. 

Note that the charging and discharging times of the proposed capacitor should be much 
smaller than the corresponding times of batteries and even smaller than times of electrolytic 
capacitors. The reason is that both batteries and electrolytic capacitors operate with the 
ion currents. In the proposed capacitor only the time of polarization of the dielectric is 
relevant. 

One can think that due to the metal infiltration the leakage current in the proposed 
capacitor should be larger than in a regular dielectric based capacitor. In principle, this 
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is true, but the effect is not strong. The conductivity of the infiltrated system near the 
percolation threshold is given by Eq.(j4]). It follows from Eq.(j5]) that s = 0.72. Then 

„t„ \ a 72 n. 28 

a[x c ) = o D a M . (26) 

Thus, the influence of metallic particles is not strong. This might be the reason why Pechar- 
roman and Moya 13 claimed a low leakage current in their experiment. 

In summary, as follows from my estimates, dielectric media with metallic nanoparticles 
that are near the percolation threshold might be promising materials for capacitors with 
large volumetric capacitance and large stored energy density. The material parameters I 
use for the estimates are far from exceptional. They do not correspond to a single specific 
material. I think that success or failure of the whole idea depends on the proper choice 
of the dielectric material. The most important requirements are high dielectric constant 
and possibility of infiltrating by metallic nanoparticles using an advanced nanotechnology. 
The requirement for the gap E g depends on whether one needs high voltage or low voltage 
capacitor. If volumetric capacitance is large enough, the stored energy density might be 
large in both cases. 

I am grateful to C. Boehme, M. Reznikov, B. I. Shklovskii, M. S. Shur, and John Worlock 
for valuable discussions and advises. 
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